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The present paper is a sequel to the paper by Karchev (2008 J. Phys.iCondens. Matter 20 325219). 
A two-sublattice ferrimagnet, with spin-si operators Sn at the sublattice A site and spin-S2 oper- 
ators S2i at the sublattice B site, is considered. Renormalized spin-wave theory, which accounts 
for the magnon-magnon interaction, and its extension are developed to describe the two ferrimag- 
netic phases (0, T*) and {T* ,Tn) in the system, and to calculate the magnetization as a function of 
temperature. 

The influence of the parameters in the theory on the characteristic temperatures Tjv and T* is 
studied. It is shown that, increasing the inter-sublattice exchange interaction, the ratio Tjv/T* > 1 
decreases approaching one, and above some critical value of the exchange constant there is only 
one phase Tjv = T*, and the magnetization-temperature curve has the typical Curie- Weiss profile. 
When the intra-exchange constant of sublattice with stronger intra-exchange interaction increases 
the Neel temperature increases while T* remains unchanged. Finally, when the magnetic order of 
the sublattice with smaller magnetic order decreases, T* decreases. The theoretical predictions are 
utilize to interpret the experimentally measured magnetization-temperature curves. 
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I. INTRODUCTION 

The present paper is a sequel to the paper [l| . A two- 
sublattice ferrimagnet, with spin-si operators Sn at the 
sublattice A site and spin-S2 operators 821 at the sub- 
lattice B site. The true magnons of a two-spin system 
are transversal fluctuations of the total magnetization 
which includes both the magnetization of the sublattice 
A and B spins. The magnon excitation is a complicate 
mixture of the transversal fluctuations of the sublattice 
A and B spins. As a result the magnons' fluctuations 
suppress, in different way, the magnetic orders on the 
different sublattices and one obtains two phases. At low 
temperature (0,T*) the magnetic orders of the A and 
B spins contribute to the magnetization of the system, 
while at the high temperature {T* ,Tm) the magnetiza- 
tion of the spins with a weaker intra-sublattice exchange 
is suppressed by magnon fluctuations, and only the spins 
with the stronger intra-sublattice exchange have non-zero 
spontaneous magnetization. 

Renormalized spin- wave theory, which accounts for the 
magnon-magnon interaction, and its extension are de- 
veloped to describe the two ferrimagnetic phases in the 
system and to calculate the magnetization as a function 
of temperature. It is impossible to require the theoret- 
ically calculated Neel temperature and magnetization- 
temperature curves to be in exact accordance with exper- 
imental results. The models are idealized, and they do 
not consider many important effects: phonon modes, sev- 
eral types of disorder, Coulomb interaction, etc. Because 
of this it is important to formulate theoretical criteria for 
adequacy of the method of calculation. In my opinion 
the calculations should be in accordance with Mermin- 
Wagner theorem 0. It claims that in two dimensions 
there is not spontaneous magnetization at non-zero tem- 
perature. Hence, the critical temperature should be equal 
to zero. It is well known that the Monte Carlo method 



of calculation does not satisfy this criteria, and "weak 
z-coupling" 3D system is used to mimic a 2D layer. It is 
difflcult within Dynamical Mean-Field Theory (DMFT) 
to make a difference between two dimensional and three 
dimensional systems. DMFT is a good approximation 
when the dimensionality goes to infinity. The present 
methods of calculation, being approximate, capture the 
basic physical features and satisfy the Mermin- Wagner 
theorem. 

There is an important difference between Neel theory 
Q and the results in the present paper. NeeVs calcula- 
tions predict a temperature T/v at which both the sublat- 
tice A and B magnetizations become equal to zero and 
T* is a temperature at which the magnetic moment has 
a maximum. 

The influence of the parameters in the theory on the 
characteristic temperatures T/v and T* is studied. It 
is shown that, increasing the inter-sublattice exchange 
interaction, the ratio Tj^/T* > 1 decreases approach- 
ing one, and above some critical value of the exchange 
constant there is only one phase T/v = T* , and the 
magnetization-temperature curve has the typical Curie- 
Weiss profile. When the intra-exchange constant of the 
sublattice with stronger intra-exchange interaction in- 
creases the Neel temperature increases while T* remains 
unchanged. Finally, when the magnetic order of the 
sublattice with smaller magnetic order decreases, T* de- 
creases. 

To compare the theoretical results and the experimen- 
tal magnetization-temperature curves one has, first of all, 
to interpret adequately the measurements. The magnetic 
moments in some materials are close to " spin only" value 
2^bS and the sublattice spins si and S2 can be obtained 
from the experimental curves. As an example I consider 
the sulpho-spinel A/nCr2iS'4_a;iS'ea; Q. On the other hand 
there are ferrimagnets with strong spin-orbital interac- 
tion. It is convenient, in that case, to consider jj coupling 
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with J"^ = L"^ + S-^ and J-^ = + S^. As an example 
I consider the vanadium spinel MnV204[l, H, 0, Hi- 

The paper is organized as follows. In Sec. II the 
model is presented and a renormalized spin-wave theory 
is worked out to calculate the magnetization-temperature 
curves for different parameters of the model. The influ- 
ence of the theory parameters on the Neel and T* tem- 
peratures is studied in Sec. III. I consider three cases: 
i) when the inter-sublattice exchange constant increases 
and all the other parameters are fixed, ii) one of the intra- 
sublattice parameters is changed and iii) when one of the 
spins decreases. Applications and analyzes of experimen- 
tal magnetization-temperature curves are given in Sec. 
IV. A summary in Sec. V concludes the paper. 



II. SPIN- WAVE THEORY 

A. Renormalized spin-wave (RSW) theory 

The Hamiltonian of the system is 

H = — Ji Sii ■ Sij — J2 • 



(1) 



where the sums are over all sites of a three-dimensional 
cubic lattice: denotes the sum over the nearest 

neighbors, <SC i,j denotes the sum over the sites of 
the A sublattice, ^ i,j 3>_b denotes the sum over the 
sites of the B sublattice. The first two terms describe 
the ferromagnetic Hciscnbcrg intra-sublatticc exchange 
Ji > 0, J2 > 0, while the third term describes the inter- 
sublattice exchange which is antiferromagnetic J > 0. To 
study a theory with the Hamiltonian Eq.([T|) it is conve- 
nient to introduce Holstein-Primakoff representation for 
the spin operators 



5y = - iSf^ = a+ .\/2si - aja^ (2) 



S 



ij - -"^i " "-j ^3 



when the sites j are from sublattice A and 

S+ = + zSi^ = -b+ ^2s2 - bp, 
S2,=Si,-zSl = -^2s2-b+h, b, (3) 
4- = -^2 + bp, 

when the sites j are from sublattice B. The operators 
a^, Qj and bj' , bj satisfy the Bose commutation relations. 
In terms of the Bose operators and keeping only the 
quadratic and quartic terms, the effective Hamiltonian 
Eq.([T]) adopts the form 



H = Ho -\- Hi 



(4) 



where 



Ho = 



siJi ^ {o.pi + apj - api - 
S2J2 i^^^^ + ^t^i - ^t^^ - ^tbj) (5) 



J ^ [sibpj + S2api - a/siS2 {afbj' + aibj)] 



+ 1-^2 E [btbph,-b,f + ibt-bp\b,] 



^ {a,bp,b, + appp,^ 



(6) 



+ \j— {afaittibj + apfa.bp - iapibpj 



and the terms without operators are dropped. 

The next step is to represent the Hamiltonian in the 
Hartrce-Fock approximation 



H w Hhf — Hci + H, 



(7) 



where 



i?c/ = 12NJisl{ui - if + 12NJ2sl{u2 - if 

+ 6NJsiS2{u-lf, (8) 



H2 = siJiUi ^ {o-tci't + apj - api - apj) 
+ S2J2U2 Y (^^^* + ^J^:> - ^J^^ - ^^^j) 



+ JuY, [^ibpj + S2api - y/siS2 {afb^ + Uibj)] 
(u) 

and N = Na = Nb is the number of sites on a sublattice. 
Equation ^ shows that the Hartrce-Fock parameters 
ui, U2 and u renormalize the intra-cxchangc constants 
Ji, J2 and the inter-exchange constant J, respectively. 

It is convenient to rewrite the Hamiltonian in momen- 
tum space representation 

Hq^Y. [^fc "fc + bpk - 7fe {apt + bkak) ] , 
fees,- 

(10) 

where the wave vector k runs over the reduced first Bril- 
louin zone Br of a cubic lattice. The dispersions are given 
by equalities 



e^. = 4si Ji itiEfe -t- 6s2 Ju 
Efc — 4s2 ^2 U2efe + 6si Ju 
7fc = 2Juy/siS2 (cos kx + cos ky + cos kz) 



(11) 
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with 

Ek = 6 - cos{kx + ky) - cos{kx - ky) - cos(fcx + kz) 
— cos{kx — kz) — cos{ky + kz) — cos(fey — kz). (12) 

To diagonalize the Hamiltonian one introduces new Bose 
fields ak, a^, Pk, (i^ by means of the transformation 



Ofc = Uk au + Vk (31 



Uk at + Vk Pk 



(13) 



hk Uk Pk + Vk a'l hi = UkPl + Vk ak 

where the coefficients of the transformation Uk and Vk 
arc real function of the wave vector k 



Uk 



(14) 



Vk = sign{"/k) „ ■ ^ 

The transformed Hamiltonian adopts the form 

[E^a+ak + E^^p+Pk + E^j , 

fee -Br 



(15) 



with new dispersions 



E2 = \ 



^{el + el)^ - 47^ - el + el 



(16) 



(4 + 4? 



ei. 



and vacuum energy 



E2^l 



+ ei)^ - 47, 
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(17) 



For positive values of the Hartree-Fock parameters and 
all values of fc G i?^, the dispersions are nonnegative 
Ek > 0, E^ > 0. For definiteness I choose Si > S2- 
With these parameters, the ak boson is the long-range 
(magnon) excitation in the two-spin system with oc 
pk^, near the zero wavevcctor, while the /3k boson is a 
gapped excitation. 

To obtain the system of equations for the Hartree-Fock 
parameters we consider the free energy of a system with 
Hamiltonian Hhf equations ^ and (fT5|) 



12Jisl{ui - 1)2 

6JsiS2{u 



1)2 



12J2SI{U2 - 1)2 



N 



(18) 



1 



E 

fceSr 



In 1 



fees. 



In 1 



where (3 = l/T is the inverse temperature. Then the 
three equations 



dT/dui = 0, dT/du2 = 0, dT/du = 
adopt the form (see the appendix) 



(19) 



Ui 



U2 



U = i - 



1 1 

6^1 N 

1 1 
6^ TV 



E 

fee Br 

E 

fee Br 



£fc 



^ifc "fc + ^'fe + «fe 



«fc< + ^fe"f + «fe 



1-- F f— f"fc< 



" + vln", + vl 



(20) 



2< + ulnl + vl 



2 / ^ ^\ (cos/c^ + cosfey + cosfc^) 



v/(£fe +4)' - 47 



where nj? and are the Bose functions of a and (3 
excitations. The Hartree-Fock parameters, the solution 
of the system of equations ([20|) . are positive functions 
of T/J, ui{T/J) > 0, U2{T/J) > and u{T/J) > 0. 
Utilizing these functions, one can calculate the sponta- 
neous magnetization of the system, which is a sum of 
the spontaneous magnetization on the two sublattices 



where 



= < Si, 



> j is from sublattice A 



(21) 



< 5*2^ > j is from sublattice B 



In terms of the Bose functions of the a and [3 excitations 
they adopt the form 



Sl 



S2 



keBr 



utnl + vlnl + 



vl 



(22) 



Uk nk + Ufc nl + w, 



/3 



feeBr 



The magnon excitation-afe in the effective theory equa- 
tion (|15p - is a complicated mixture of the transversal fluc- 
tuations of the A and B spins. As a result the magnons' 
fluctuations suppress in a different way the magnetization 
on sublattices A and B. Quantitatively this depends on 
the coefficients Uk and Vk in equations (j22p . At char- 
acteristic temperature T* spontaneous magnetization on 
sublattice B becomes equal to zero, while spontaneous 
magnetization on sublattice B is still nonzero. Above 
T* the system of equations ([20)1 has no solution and one 
has to modify the spin-wave theory. The magnetization 
depends on the dimensionless temperature T/J and di- 
mensionless parameters si, S2^ Ji/J and J2I J ■ For pa- 
rameters Sl = 1.5, S2 ^ \t J\l J = 0.94 and J2/J = 0.01 
the functions M^^TjJ) and {T j J) are depicted in 
figure 1. The upper (blue) line is the sublattice A mag- 
netization, the bottom (red) line is the sublattice B mag- 
netization. 
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FIG. 1: (color online) The spontaneous magnetization M^- 
upper (blue) line and Af^-bottom (red) line as a function of 
T/J for parameters si = 1.5, S2 — 1, Ji/J — 0.94 and 
J2/J = 0.01. T* is the temperature at which sublattice B 
magnetization becomes equal to zero 



B. Modified RSW theory 

Once suppressed, the sublattice B magnetization can- 
not be restored increasing the temperature above T*. To 
formulate this mathematically we modify the spin-wave 
theory using the idea of a description of the paramag- 
netic phase of 2D ferromagnets {T > 0) by means of 
modified spin- wave theory jl^, and its generalization 
We consider a two-sublattice system and, to enforce 
the magnetization on the two sublattices to be equal to 
zero in paramagnetic phase, we introduce two parame- 
ters Xa and As The new Hamiltonian is obtained 
from the old one equation ([T|) by adding two new terms: 



(23) 



In momentum space the new Hamiltonian adopts the 
form 



keBr 



where the new dispersions are 



(24) 



il=el + X,, + (25) 



Utilizing the same transformation equations (jl3p with 
parameters 



el 



(26) 



Vk = sign{jk) 



one obtains the Hamiltonian in diagonal forma 



(27) 



k£Br 



where 



El 



{el + 4)2 - 47, 



+ £ 



(^1 + 4) 



b\2 



(28) 



{el + eiY - 472 



It is convenient to represent the parameters Ai and A2 
in the form 

Ai ~ 6Jus2{ni — 1), A2 = 6Jusi{^2 — !)■ (29) 



In terms of the new parameters /ii and /i2 the dispersions 
3 form 

AsiJiUi ek + 6s2 Jw/ii 



el and e^. adopt the form 



ei. 



el 



(30) 



4S2 J2 U2 £fc + GsiJufl2 



They are positive {el > 0, £^ > 0) for all values of the 
wavevector k, if the parameters fii and fj,2 are positive 
(/ii > 0, /i2 > 0). The dispersions Eg. ipS)) are well de- 
fined if square-roots in equations ([28]) are well defined. 
This is true if 



fllfl2 > 1. 



(31) 



The Pk excitation is gapped {E^ > 0) for all values of 
parameters /xi and ^2 which satisfy equation (pij) . The 
a excitation is gapped if /ii/i2 > 1, but in the particular 
case 



^ = II [^fe + ei b+bk - Ik {bkak + b+al)] 



Q — 0, and near the zero wavevector 

Et « Pk^ 
with spin-stiffness constant 



8(s2-^2M2Ail + s\JiUi^2) + 2siS2Ju 
(Sl^2 - S2M1) 



(32) 



(33) 



(34) 



In the particular case equation (|32p boson is the long- 
range excitation (magnon) in the system. 

We introduced the parameters Ai and A2 (/ii , /i2) to en- 
force the sublattice A and B spontaneous magnetizations 
to be equal to zero in the paramagnetic phase. We find 
out the parameters /ii and /i2, as well as the Hartree-Fock 
parameters, as functions of temperature, solving the sys- 
tem of five equations, equations (j20p and the equations 
= = 0, where the spontaneous magnetizations 
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have the same representation as equations ((22|) but with 
coefficients iik, Vk: and dispersions E'j^, in the ex- 
pressions for the Bose functions. The numerical calcula- 
tions show that for high enough temperature ^.1^2 > 1- 
When the temperature decreases the product ^1^2 de- 
creases, remaining larger than one. The temperature at 
which the product becomes equal to one (/ii/i2 = 1) is 
the Neel temperature. Below T/v, the spectrum contains 
long-range (magnon) excitations, thereupon /ii/i2 = 1- It 
is convenient to represent the parameters in the following 
way: 

Ml = Mi M2 = (35) 

In the ordered phase magnon excitations are the ori- 
gin of the suppression of the magnetization. Near the 
zero temperature their contribution is small and at zero 
temperature sublattice A and B spontaneous magneti- 
zation reach their saturation. On increasing the tem- 
perature magnon fluctuations suppress the sublattice A 
magnetization and sublattice B magnetization in differ- 
ent ways. At T* the sublattice B spontaneous magne- 
tization becomes equal to zero. Increasing the tempera- 
ture above T* , the sublattice B magnetization should be 
zero. This is why we impose the condition M^{T) = 
if T > T*. For temperatures above T* , the parameter /i 
and the Hartree-Fock parameters are solution of a sys- 
tem of four equations, equations ([^D|) and the equation 
= 0. The Hartree-Fock parameters, as a functions of 
temperature T/ J, are depicted in figure 2 for parameters 
si = 1.5, S2 = 1, Ji/J = 0.94 and J2/J = 0.01. The 
vertical dotted (green) line corresponds to T* / J. 
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FIG. 3: (color online) fi{T/J) for parameters si = 1.5, S2 = 
1, Ji/J = 0.94 and J2/J = 0.01. The vertical dotted (green) 
line corresponds to T* / J , while (red) dashed lines to Tn / J 
and fi{TN/J). 



hmction of the temperature. Above T*, the magneti- 
zation of the system is equal to the sublattice A magne- 
tization. For the same parameters as above the functions 
A/^(r/J) and M^{T/J) are depicted in figure 4a. The 
upper (blue) line is the sublattice A magnetization, the 
bottom (red) line is the sublattice B magnetization. The 
total magnetization M = + is depicted in fig- 
ure 4b. 
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FIG. 2: (color online) Hartree-Fock parameters u\,U2 and u 
as a function of T/J for si = 1.5, S2 = 1, Ji/J = 0.94 and 
J2/J = 0.01. The vertical dotted (green) line corresponds to 
TVJ 

The function fi{T/J) is depicted in figure 3 for the 
same parameters. 

We utilize the obtained function fi{T), Ui{T), U2{T), 
u(T) to calculate the spontaneous magnetization as a 



III. Tn and T* dependence ON MODEL'S 
PARAMETERS 



The existence of two ferromagnetic phases (0, T*) and 
(T*,r/v) is a generic feature of two spin systems. The 
characteristic temperatures Tn and T* strongly depend 
on the parameters of the model. Intuitively, it is clear 
that, if the inter-exchange is much stronger than intra- 
exchanges, the ferromagnetic order sets in simultane- 
ously on both sublattices. This is not true, if inter- 
exchange is not so strong. To demonstrate this I study 
a system with sublattice A spin si = 1.5, and sublat- 
tice B spin S2 = 1- For parameters Ji/J = 0.5 and 
J2/J = 0.005 the magnetization-temperature curve is 
depicted in FIG. 5 curve "c". The ratio of the charac- 
teristic temperatures equals T/v/T* = 1.722. Increasing 
the inter-exchange coupling, Ji/J = 0.3, J2/J = 0.003 
(curve "b"), the ratio decreases T^/T* = 1.229, and 
above some critical value of the inter-exchange constant 
Ji/J = 0.05, J2/J = 0.0005 NeeVs temperature becomes 
equal to T*. There is only one ferromagnetic phase, and 
magnetization-temperature curve "a" is a typical Curie- 
Weiss curve. Despite this the system does not describe 
fcrromagnet, because the spin wave excitations are su- 
perposition of the sublattice A and B spin excitations. 
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FIG. 5: (color online) The magnetization 2 M'^ + 2 Af^ as 
a function of T/J for si = 1.5 and S2 = 1, curve a:Ji/J = 



0.05, J2/J = 0.0005, curve b:Ji/J 
curve c:Ji/J = 0.5, J2/J = 0.005. 



= 0.3, J2/J = 0.003, 



FIG. 4: (color online) a) The sublattice A spontaneous mag- 
netization A/'^-upper (blue) line and sublattice B sponta- 
neous magnetization A'/^-bottom (red) line as a function of 
T/J for parameters si = 1.5, S2 = 1, Ji/J = 0.94 and 
J2/J = 0.01. 

b) The total spontaneous magnetization + . T* /J- 
vertical dotted (green) line 



Next, I consider a system with sublattice A spin si = 
1.5, and sublattice B spin S2 = 1 . The ratio of sublat- 
tice B exchange constant J2 and inter-exchange constant 
J is fixed j2 = J2I J = 0.01, while the ratio ji = Ji/J 
varies. When the sublattice A exchange constant Ji in- 
creases ji = Ji/J = 0.64, 0.84, 0.94, the magnetization- 
temperature curve at temperatures below T* does not 
change. There is no visible difference between T* tem- 
peratures for the three values of the parameter Ji/J. 
The difference appears when the temperature is above 
T* . Increasing sublattice A exchange constat increases 
the Neel temperature. The three curves are depicted in 
figure 6. 

Finally, I consider three systems with equal exchange 
constants Ji/J = 0.4, J2/J = 0.004 and sublattice A 
spin si = 4, but with three different sublattice B spins 
(figure 7). The calculations show that decreasing the 
sublattice B spin decreases T* temperature, increases the 
maximum of magnetization at T* and zero temperature 
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FIG. 6: (color online) The magnetization 2 + 2 M'^ as a 
function of T/J for si = 1.5, S2 = 1, 32 ~ Ji/J ~ 0.01 and 
three values of the parameter ji = Ji/J; ji = 0.94 (black) 
squares, ji = 0.84 (red) circles, ji — 0.64 (blue) triangles 



magnetization. 
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c (blue). The temperature and magnetization axis are 
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FIG. 7: (color online) The magnetization 2M^ + 2M^ as 
function of T/J for Ji/J = 0.4, J2/J = 0.004, si = 4, and 
S2 = 0.5-curve a (green), S2 = 1.5-curve b(red), S2 = 2.5- 
curve c (black) 



IV. THEORY AND EXPERIMENT 
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FIG. 8: (color online) The magnetization 2M^ + 2M^ as 
function of T/J for Ji/J = 0.47, J2/J = 0.001, si = 1.5, 
and S2 = 1-curve a (black), S2 = 0.7-curve b(red), S2 = 0.4- 
curve c (blue) 



A. Sulpho-spinel MnCr2Si-xSex 

The sulpho-spinel MnCr2Si-xSex has been investi- 
gated by measurements of the magnetization at 15.3fcOe 
as a function of temperature (figure 94 in j3l). The max- 
imum in the magnetization versus temperature curve, 
which is typical of AInCr2Si {x = 0), increases when 
X increase, and disappears at x = 0.5. The Neel tem- 
perature decreases from 74 A' at a; = to 56 A' aX x — 2. 
The authors' conclusion is that the observed change of 
the magnetic properties is attributed to a decrease of the 
strength of the negative Mn?^ — Cr^^ superexchange 
interaction with increasing Se concentration. 

We obtained, see figure 5, that the maximum of the 
magnetization is at T* . Above T* the magnetization of 
the system is equal to the magnetization of sublattice A 
spins. If we extrapolate this curve below T* down to zero 
temperature we will obtain a value close to 2sifiB, where 
si is the spin of the sublattice A spin operators. The 
experimental figures Q show that extrapolations give one 
and the same result for all values of x. One can accept 
the fact that the 5*6 concentration do not influence over 
the value of sublattice A spin and si = 1.5. 

Below T* the magnetization is a sum of sublattice 
A and B magnetization. Hence, the magnetization at 
zero temperature is equal to 2(si — S2)a*s- Therefore, 
one can determine the sublattice B spin S2- The re- 
sults of the theoretical calculations of magnetization, in 
Bohr magnetons, are depicted in figure 8 for parameters 
si = 1.5, Ji/J = 0.47, J2/J = 0.001 and 53 = 1-curve 
a(black); S2 — 0.7-curve b (red), and S2 = 0.4-curve 



chosen in accordance with experimental figure. Compar- 
ing figure 94 in [3| and figure 7 in the present paper, one 
concludes that the effective sublattice B spin S2 decreases 
with increasing Se concentration, and this is the origin 
of the anomalous temperature variation of magnetiza- 
tion. The figure 8 shows that the present calculations 
capture the essential features of the system; increasing 
the Se concentration (decreasing S2) leads to a decrease 
of Neel temperature, T* temperature decreases too, and 
the maximum of the magnetization increases. Compar- 
ing the figure 8 in the present paper and figure 5 in [l[ 
one realizes the importance of the present method of cal- 
culation for adequate reproducing the characteristic tem- 
peratures T/v, T*, and the shape of the magnetization- 
temperature curves. 



B. Vanadium spinel MnV204, 

The spinel MnV204 is a two-sublattice ferrimagnet, 
with site A occupied by the Mn?^ ion, which is in the 
high-spin configuration with quenched orbital angular 
momentum, which can be regarded as a simple s = 5/2 
spin. The B site is occupied by the V^^ ion, which takes 
the 3c?^ high-spin configuration in the triply degenerate 
t2g orbital and has orbital degrees of freedom. The mea- 
surements show that the setting in of the magnetic order 
is at Neel temperature Tjv = 56. 5A' 0] and that the 
magnetization has a maximum near T* = 53.5Ar. Be- 
low this temperature the magnetization sharply decreases 
and goes to zero when temperature approaches zero. 
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We consider a system which obtains its magnetic prop- 
erties from Mn and V magnetic moments. Because of the 
strong spin-orbital interaction it is convenient to consider 
jj coupling with = and = -t- S^. The sub- 
lattice A total angular momentum is Ja = sa = 5/2, 
while the sublattice B total angular momentum is = 
Ib + sb, with Ib = 5, and ss = 1 Then the g- factor 
for the sublattice A is gA = "2, and the atomic value of 
the gB is Ob ~ J- The sublattice A magnetic order is 
antiparallel to the sublattice B one and the saturated 
magnetization is tr = 2| — |4 = 0, in agreement with 
the experimental finding that the magnetization goes to 
zero when the temperature approaches zero. The Hamil- 
tonian of the system is 



H ^ - 



A -tA 



KB ^ J 



B tB 



(36) 



The first two terms describe the ferromagnetic Heisen- 
berg intra-sublattice exchange ka > 0,kb > 0, while 
the third term describes the inter-sublattice exchange 
which is antifcrromagnetic k > 0. To proceed we use 
the Holstein-Primakoff representation of the total angu- 
lar momentum vectors Jj^(a^, a^) and J^(6^, bj), where 
a'j', Qj and b'j' , bj are Bose fields, and repeat the calcula- 
tions from sections II and III. The magnetization of the 
system gA + gB as a function of the tempera- 
ture is depicted in figure 9 for parameters ka/k = 0.45 
and kb/k = 0.001. The parameters are chosen so that 
the calculations to reproduce the experimental value of 
the ratio Tn/T*. 
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FIG. 9: (color online) The magnetization qa + Qb as 
a function of T/k for parameters ka/k = 0.45 and kb/k = 
0.001. 

The profile of the magnetization-temperature curve is 
in a very good agrcment with the experimental zero- 
field cooling (ZFC) magnetization curves 0, 0]. The 



anomalous temperature dependence of the magnetization 
is reproduced, but there is an important difference be- 
tween the interpretation of the experimental results in 
@, H, 0j H, 0, and the present theoretical results. In 
the experimental papers Tjv is the temperature at which 
both the Mn and V magnetization become equal to zero. 
The present theory predicts two phases: at low temper- 
atures (0, r*) sublattice Mn magnetization and sublat- 
tice V magnetization contribute to the magnetization of 
the system, while at high temperatures {T*,Tn) only 
Mn ions have non-zero spontaneous magnetization. The 
vanadium sublattice magnetization set in at T*, and ev- 
idence for this is the abrupt decrease of magnetization 
below T*, which also indicates that the magnetic order 
of vanadium electrons is anti-parallel to the order of Mn 
electrons. 

For samples cooled in a field (PC magnetization) the 
field leads to formation of a single domain and, in addi- 
tion, increases the chaotic order of the spontaneous mag- 
netization of the vanadium sublattice, which is antipar- 
allel to it. As a result the average value of the vana- 
dium magnetic order decreases and does not compensate 
the Mn magnetic order. The magnetization curves de- 
pend on the applied field, and do not go to zero. For 
a larger field the (FC) curve increases when tempera- 
ture decreases below Neel temperature . It has a max- 
imum at the same temperature T* < Tjv as the ZFC 
magnetization, and a minimum at Tj* < T*. Below Tj* 
the magnetization increases monotonically when temper- 
ature approaches zero. 

The experiments with samples cooled in field (FC mag- 
netization) provide a new opportunity to clarify the mag- 
netism of the manganese vanadium oxide spinel. The 
applied field is antiparallel with vanadium magnetic mo- 
ment and strongly effect it. On the other hand, the ex- 
periments show that there is no difference between ZFC 
and FC magnetization curves when the temperature runs 
over the interval (T*,Tjv) ^M- They be gin to diverge 
when the temperature is below T*. This is in accor- 
dance with the theoretical prediction that the vanadium 
magnetic moment does not contribute the magnetization 
when T > T* and T* is the temperature at which the 
vanadium ions start to contribute the magnetization of 
the system. Because of the strong field, the two vana- 
dium bands are split and the magnetic moment of one of 
the t2g electrons is reoriented to be parallel with the field 
and magnetic order of the Mn electrons. The description 
of this case is more complicate and requires three mag- 
netic orders to be involved. When T* < T < T/v only 
Mn ions have non zero spontaneous magnetization. At 
T* vanadium magnetic order antiparallel to the magnetic 
order of Mn sets in and partially compensates it. Below 
T* the reoriented electron gives contribution, which ex- 
plains the increasing of the magnetization of the system 
when the temperature approaches zero. A series of ex- 
periments with different applied field could be decisive 
for the confirmation or rejection of the T* transition. In- 
creasing the applied field one expects increasing of T* 
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and when the field is strong enough, so that all vana- 
dium electrons are reoriented, an anomalous increasing 
of magnetization below T* would be obtained as within 
the ferromagnetic phase of UGe2 • 

V. SUMMARY 

In summary, I have worked out a rcnormalizcd spin- 
wave theory and its extension to describe the two phases 
(0,r*) and (T*,Tn) of a two sublatticc fcrrimagnet. 
Comparing the figure 4 in the present paper and figure 4 
in [l| and figure 8 in the present paper and figure 5 in [l[ 
one becomes aware of the relevance of the present calcula- 
tions for the accurate reproduction of the basic features 
of the system near the characteristic temperatures Tn 
and T*. 

The present theory of ferrimagnetism permits to con- 
sider more complicate systems such as CeCrSb^ com- 
pound or the spinel Fe304 which are two sublattice 
ferrimagnets but with three spins. 
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APPENDIX A 

To make more transparent the derivation of the equa- 
tions for the Hartree-Fock parameters Eg. ipO)) I consider 
the first term (the sublattice A term) in the Hamiltonian 
of the magnon-magnon interaction Eq.®. To write this 
term in the Hartree-Fock approximation one represents 
the product of two Bose operators in the form 

o-t"-j = "i^flj - < > + < a+a-i > (Al) 

and neglects all terms {a^aj — < afaj in the four 
magnon interaction Hamiltonian. The result is 

ia+flja+a, w - < a+a^ >< a+a, > 

+ < d^aj > afai + a'^aj < a'^Ui > 

-a^aiO^aj w - < a+Oj >< a^aj > 

+ < a+Ci > a^aj + ai < a'j'aj > 

-a^ajafaj w — < a+a,- >< a^a_, > (A2) 
2 

+ < ajuj > a^aj + a^aj < a^aj > 



+ < afai > a^aj + a'^at < d^a^ > 

— < a^aj >< a^Oi > 

-I- < alaj > a^ai + a^Ui < afuj > 
The Hartree-Fock approximation of the sublattice A part 
of the Hamiltonian of magnon-magnon interaction reads 

i Ji [ata+iai - aj f + («+ - ap^maj] 

« l2NJisl (mi - if (A3) 
+ JiSi {ui — 1) ^ (af'^i + — a^ai — afuj) 

where the Hartree-Fock parameter ui is defined by the 
equation 

"1 = 1- E < ^t^k > (A4) 

Combining the sublattice A part of the Hamiltonian 
Eq.([5]) (the first term) and Ea. (|A3p one obtaines the 
Hartree-Fock approximation for the sublattice A part of 
the Hamiltonian 

« l2NJisl {ui - if (A5) 
-I- Jisiui E {o'td + aljUj — a^Ui — afa^) 

In the same way one obtains the Hartree-Fock approx- 
imation of the sublattice B and inter sublattices parts 
of the Hamiltonian. The result is the Hhf Hamiltonian 
Eas. (1715151) . 

To calculate the thermal average < a'^a^ >, in the 
Ea. (jA4[) . one utilizes the Hamiltonian Hhf- Therefor, 
the matrix element depends on the Hartree-Fock param- 
eters, and equation (|A4p is one of the self consistent equa- 
tions for these parameters. 

The matrix element can be represented in terms of 
ak{a+) and /?fe(/?+) Eq.^ 

< a+Uk > = ui nl + vl n'l + vl (A6) 

where =< a^Uk >, "■f =< Pt > are the Bose 
functions of a and (3 excitations. Substituting the ther- 
mal average in Ea. (|A4p with Ea. (|A6p . one obtains that 
equation (|A4p is exactly the first equation of the system 
Ea. ipO]) which in turn is obtained from the first of the 
equations 
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